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ABSTRACT 

The site and mechanism of gamma-ray burst (GRB) prompt emission is still un- 
known. Although internal shocks have been widely discussed as the emission site 
of GRBs, evidence supporting other emission sites, including the closer-in photo- 
sphere where the fireball becomes transparent and further-out radii near the fireball 
deceleration radius where magnetic dissipation may be important, have been also 
suggested recently. With the successful operation of the GLAST experiment, prompt 
high energy emission spectra from many GRBs would be detected in the near future. 
We suggest that the cut-off energy of the prompt emission spectrum from a GRB de- 
pends on both the fireball bulk Lorentz factor and the unknown emission radius from 
the central engine. If the bulk Lorentz factor could be independently measured (e.g. 
from early afterglow observations), the observed spectral cutoff energy can be used to 
diagnose the emission site of gamma-rays. This would provide valuable information 
to understand the physical origin of GRB promp emission. 



1 INTRODUCTION 

The em ission of photons in the prompt phase of a GRB may last from less than a second to hundreds of seconds (see e.g. lMeszarosI 
200d) for a recent general re view on GRBs). The exac t physical me chanism and the emission site of the observed prompt GRB 
emissions are still unknown fchang & Meszaros! |2004| ; lzhang||2007h . The physical processes which may lead to this emission 
include synchrotron/jitter emission, inverse Compton scattering or a combination of thermal and non-thermal emission components. 
Internal shocks have been widely discussed in the literature as the possible emission site of GRB prompt emission (Rees & Meszaros 
1994; Meszaros et al. 1994; Kobayashi et al. 1997; Daigne & Mochkovitch 1998; Pilla & Loeb 1998; Panaitescu & Meszaros 2000; 
Lloyd & Petrosian 2000; Zhang & Meszaros 2002a; Dai & Lu 2002; Pe'er & Waxman 2004, 2005; Pe'er et al. 2005; Gupta & Zhang 
2007). Within this model the radius of emission from the central engine (r) is related to the bulk Lorentz factor (T) and the variability 
time (t„) through r = Y 2 ct v . For typically observed values of bulk Lorentz factor V ~ 300 and variability time t v ~ 0.01 sec the 
emission radius is r ~ 3 x 10 13 cm. However, the internal shock origin of GRB prompt emission is not conclusive. The baryonic 
or pair photospheres of the GRB fireball have been argued to be another possible emission site of prompt GRB emission (e.g. 
Thompson 1994; Meszaros & Rees 2000; Meszaros et al. 2002; Ryde 2005). This model has the merit of potentially reproducing 
some observed empirical correlations among GRB prompt emission properties (e.g. Rees & Meszaros 2005; Ryde et al. 2006; 
Thompson et al. 2007; cf. Zhang & Meszaros 2002). On the other hand, an analysis of Swift early afterglow data led Kumar et 
al. (2007) to conclude that the prompt emission site is between 10 15 -10 16 cm from the central engine. This emission site is too 
large for typical internal shocks but too small for external shockfl Emission at this radius may be related to magnetic dissipation 
(e.g. Spruit et al. 2001; Drenkhahn & Spruit 2002; Zhang & Meszaros 2002; Giannios & Spruit 2007). In both of the above two 



1 This radius can be still accommodated within the internal shock picture if the typical variability time scale is a significant fraction of the burst 
duration. Liang et al (2006) discovered that if the steep decay segment observed in Swift X-ray afterglows is due to the curvature effect of the high 
latitude emission with respect to the line of sight (Kumar & Panaitescu 2000; Zhang et al. 2006), the required time zero point to usually leads the 
beginning of the steep decay (t p ), and (t p — to) (effectively the variability time scale for the internal shock scenario) is a significant fraction of the 
burst duration. The internal shock scenario therefore can be still consistent with Kumar et al.'s analysis. 
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non-internal shock models for GRB prompt emissions, it is possible to argue that the emission radius r could be in principle related 
to the Lorentz factor Y and variability t v in a non-trivial (e.g. other than r = F ct v ) manner. For example, in the photosphere 
model, r is defined by the optically-thin condition, and is not directly related to t v which is related to the time history of the GRB 
central engine. In the magnetic dissipation model, if the energy dissipation occurs locally (i.e. the emission region scale is much 
smaller than the emission radius), it is possible to have r > Y 2 ct v - In general, it is reasonable to treat r as an independent quantity 
with respect to Y and t v . 

High energy photons produced in the prompt emission region are expected to interact with lower energy photons before 
escaping as a result of two photon attenuation. In general, the internal optical depth of 77 interactions depends on F, t v , and the 
radius of the emission region. Fraditionally, internal shocks have been taken as the default model of GRB prompt emission, and 
the pair attenuation optical depth has been expressed as a function of Y and t v only (e.g. Piran 1999; Fithwick & Sari 2001). Fhe 
pair attenuation process is expected to leave a cutoff spectral feature in the prompt emission spectrum, and detecting such a spectral 
cutoff by hi gh energy missions such as GFASF ha s been discussed as an important method to estimate the bulk Forentz factor F of 
the fireball tearing & Hardindll997l ; lBaring||2006l) . Fhe issue of unknown emission radius r makes the picture more complicated. 
It is no longer straightforward to estimate Y with an observed spectral cutoff energy. On the other hand, there are other independent 
methods of estimating Y using early afterglow (e.g. Sari & Piran 1999; Zhang, Kobayashi & Meszaros 2003) or prompt emission 
(Pe'er et al. 2007) data and there have been cases of such measurements (Molinari et al. 2007; Pe'er et al. 2007). It is then possible 
to use the observed spectral cutoff energy to diagnose the unknown GRB emission site if Y is measured by other means. In this 
paper we release the internal shock assumption and re-express the cutoff energy more generally as a function of r and Y. We then 
discuss an approach of diagnosing the GRB prompt emission site using the future cutoff energy data retrieved by GFASF and 
other missions. Fately Murase & Ioka (2007, see also Murase & Nagataki 2006) also independently discussed to use the pair cutoff 
signature to diagnose whether the emission site is the pair/baryonic photosphere. We discuss this topic more generally to diagnose 
any emission site. We emphasise that our method can be used to constrain r only when a clear cut-off is observed in the high energy 
photon spectrum from a GRB. 



2 PARAMETRIZATION OF INTERNAL OPTICAL DEPTH 

The cross section of two-photon interaction can be generally expressed as (Gould & Schreder 1967) 

(1) 
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where <jt is the Thomson cross section, b = [1 — (E' th /E' yi J] 1 ' 2 is the center of mass dimensionless speed of the pair produced, 
E h , E lt and 8' are the high- and low-energy photon energies and their incident angles in the comoving frame of the GRB ejecta. 
The threshold energy of pair production for a high energy photon with energy E'^ h is 
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In the comoving frame, the relative velocity of the high energy and low energy photons along the direction of the former is c(l — 
cos 6'). For an isotropic distribution the fraction of low energy photons moving in the differential cone at an angle between 8' and 
(8' + d8') is i sin 9'd8'. The inverse of the mean free path for jh 7; interactions IZ. „ (EL ) can be calculated as 

where 
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— ,„i — is the specific number density of low energy photons of the GRB in the comoving frame. The observed low energy 
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photon spectrum (per unit energy per unit area) d ^ a gjjj from a GRB pulse can be used to estimate — dE , — , i.e. 

d nil {EL a ) _ dl dN° { (E° t ) dE° ( 



dE!y t r 2 A' dE° t dA dE'~ n 



© 2007 RAS, MNRAS 000,[T]-?? 



Diagnosing GRB Prompt Emission Site with Spectral Cut-Off Energy 3 



where A' is the comoving width of the shell at the radius r from the central engine. We notice that the expression of A' is function of 
radius (e.g. Zhang & Meszaros 2002b): A' = r for r < r c ; A' = r c = Fct v for r c < r < r s \ and A' = r/T for r > r s = T 2 ct v . 
Throughout the paper the superscript "o" denotes the quantities measured in the observer's rest frame. 
The comoving distance of the source is 

*. = wf i =, (6) 
H " Jo y/n A + n m (i + xy 

which is related to the luminosity distance through — d z (l + z), where z is the redshift of the source. Here Ho = 
71km s _1 Mpc -1 Qa = 0.73 and Q m — 0.27 are adopted in our calculations. The observed fluence is usually a broken power 
law with a break energy of the order of MeV (Band et al. 1993). We assume that the two spectral indices below and above the break 
energy are — /3i and — /?2, respectively, usually with f3\ ~ 1 and P2 > 2. We model the observed photon flux as 

dK(K) _ f of K.- 01 E° n <E° M m 

dE° t dA J \ E° lM ^~^E°-^ E° l >E° M 

The break energy in the photon spectrum in the observer's frame is denoted by E° br . We define a new variable following the 
procedure discussed in Gould and Schreder (1967) to reduce the number of integrals 
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with so — ^ c 'iyl > an d O — ^ (1 — c °s ). As 6 = (1 — 1/s) ' , the pair production cross section can be expressed as a function 
of s. It is then possible to write Eq.([3} as 



r 1 i 



where 

Q[s (S; i )] = J sa(s)ds, (10) 

and a(s) = ^ "H^"' ■ F° r moderate values of s we use er(s) ~ 1 and the expressions for <3[s (E^)] becomes (sg - l)/2. 
Substituting for Q[so(£ , 7 i )] in Eq.([9l( we derive the final expression for l^ hll (E' lh ). Finally, the internal optical depth Ti nt (E° h ) 
is the ratio of width of 7^7; interaction region and the mean free path of their interaction. In most cases, this is simply 

T int (E° h ) = A'/l^ hyi (E; h ). (11) 

In this case, comparing Eqs.((5} and i ll It suggests that the concrete expression of A' does not enter the problem since it is 
cancelled out in the expression of Tj n t( E° u ). In the photosphere models that invoke a continous wind from the central engine 
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the expression ot Ti n t\cj lh ). in tne photosphere models that invoke a 
s 2000: lGiannioi 2006). however, this expression should be modified as 



nnt{E° h )=r/[Tl^ l {E° h )). (12) 

We therefore also consider such a case. 

In the following we discuss three cases. The analytical expressions below are only valid for Eg. ill It. For the case of continuous 
photosphere models (Eg. 1121 ). the analytical expressions for r < r is are more complicated, and we only present the numerical 
results in Fig.l, where A' = Tct v is used. 

Case (I): Both the cutoff energy E° c (defined by Ti nt (E° c ) — 1 at which energy the observed spectrum significantly deviates from 
the power-law extension of the low energy spectrum) and the threshold energy E° { th (for E° c ) are above the break energy in the 
observed photon spectrum E° br , i.e. E° br < E° l th < E° c (or E^ br < E! yi th < E' lc ). The expression for the optical depth is 
the simplest for this case 

Ai(E° ) / r \ 2 " 2 ' 3 2 

n nt (^J = ^2( T ^-) (13) 
where 
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radius of prompt emission is r = F 2 ct° /(l + 2) where t° is the observed variability time scale. Substituting this expression of r 
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a™ 1 dE° dl = /i£°r Wlth f? = / #°,(,r has been assumed for E° t > E° M . In the case of internal shocks the 
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in Eq.lll3t one gets n„t(E° ) oc Y~ 2 ~ 2/32 , which is consistent with lLithwick & S aril 12001). Our expression is more generic with 
r being a free parameter. 

Case (II): If the cutoff energy is still above the break energy, but the threshold energy for pair production is below the break energy, 
i.e. (E° l th < E° br < E° c ), the expression for internal optical depth is more complicated. Making use of Eq.0, one gets 

A 2 (E° h ) / r \ 2-2/3 



n nt {E° h ) = 
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for /3i 7^ 1, where E° h br is the energy of the high energy photons that interact with the break-energy photons at the threshold 

, e 2 c 4 . Equation d 1 6fc can be reduced to Eq. J14t when f3i = j3%. For 



r 
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condition, which is defined by E° hbr E° br = 
f3\ — 1, Tint does not depend on V, and one has 
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In this case T is not needed to infer r. 

Case (III): In more extreme cases, usually with a low enough Lorentz factor, one could have the cutoff energy below the break 
energy, i.e. E° uth < E° c < E° >br . 

In this regime, we still use Eq.© to calculate the internal optical depth, but effectively one can place the upper limit of the 
integration as the break energy, since above the break energy the photon flux falls off rapidly. This gives 

A 3 (E°)/ T \ 2 — 2/3i 



Tint(E° h ) = 
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where 
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(19) 



To summarize all three cases, the radius of the prompt emission can be calculated in terms of cutoff energies by making the 
internal optical depths (Eqs J13ll5ll8l ) to unity 



+ z 



i-fj j 



(20) 



where j = 2 for i = 1, and j — 1 for i — 2,3. In practice, from the observed low energy photon spectrum it is possible to derive 
the low energy photon spectral parameters (including f}\, /3 2 , E° br , E° c , etc.) and to identify one applicable case among the three 
cases discussed. If the burst redshift z is measured from afterglow observations, and if the GRB bulk Lorentz factor is measured or 
constrained independently with other methods (e.g. Zhang et al. 2003; Molinari et al. 2007; Pe'er et al. 2007), one can estimate the 
GRB emission radius using Eq.d20l>. For the case (I) and if 02 ~ 2, a very simple expression of r is available according to Eq.<ll3b 



dr 



m„c 2 r 
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which could be used to quickly estimate r with the data. 



3 CASE STUDY AND DISCUSSION 

We have generalized the expression of the cutoff energy of prompt GRB spectrum to include the dependences of both V and r. 
We suggest that the information of this cutoff energy is useful to diagnose the unknown location of gamma-rays. We discuss three 
cases, and derive a general expression of r (Eg. 1201 ). In Fig. 1(a) we present the contours of E° o in the r — T plane (detailed 
model parameters are listed in the figure caption). It is straightforward to see that E° a carries the information of both V and r, and 
unless the internal shock model is assumed, one cannot constrain T by E° a . Figure 1(b) indicates that by knowing Y from other 
measurements, one could diagnose r with the E° a information. GLAST will be launched in early 2008, and LAT on board will be 
able to measure the pair cutoff signature for many GRBs. With the coordinated observations with Swift and ground-based follow up 
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10 2 1()3 "" 1 10 10 3 10 5 

r ' ^;(MeV) 

Figure 1. (a) Contours of E° c in the r - Y plane. From top to bottom, E° c =50GeV, 5GeV, 0.5GeV, 5()MeV, 5MeV and 0.5MeV, respectively. 
Model inputs: E° br = lMeV, observed flux is 0.6MeV _1 cm- 2 at E° br , z = 1, fix = 1.5 and /3 2 = 2.2. The solid lines are derived from 
Eq.(TTJ. The dotted lines are derived from Eq. )12t for r < ri s . The dashed line represents the internal shock r — Y relation, and the variability 
time t v = O.Olsec is assumed, (b) Inferring r from the observed cutoff energy. From top to bottom, Y = 100, 200, 400, 600, 800, and 1000, 
respectively. Other parameters and line styles are the same as (a). The crosses on the lines for Y = 400 to 1000 are for the internal shock model 
with t v = O.Olsec. 



observations (to obtain z and Y information), one may be able to more precisely diagnose the emission site of gamma-rays, which 
is so-far subject to debate. 



We notice that the maximum observed photon energy (e max ) mentioned in the paper bv lLithwick & Sari (2001) is not the 
same as the cut-off energy {EZ ) discussed in our paper. They discussed the maximum observed energy defined by the detector 
bandpass and sensitivity, and used that energy to derive the lower limit of T. Within that context, they discussed two possible 
limits defined by pair attenuation (limit A) and Compton scattering by pairs (limit B), respectively. Here we discuss the physical 
cutoff e nergy, so we can i nfer th e actual value (not lower limits) of Y and r (our new addition). By definition, by regarding E° a as 



in llithwick & Sari 200 lh . their e maXjan is just ethick, and their self-annihilation energy e sc if, an is always smaller than or at 



most equal to em a x by definition (otherwise photons cannot be attenuated at e max = E° c ). As a result, the limit B discussed by 



Lithwick & S aril a200lh is never relevant in our context. 

We have assumed that the photon field is isotropic in the comoving frame. In some models (e.g. Lyutikov & Blandford 2003; 
Thompson et al. 2007) the emitters are moving fast in the comoving frame of the bulk flow. This will introduce anisotropy of the 
photon field in the comoving frame. Suppose the relative bulk Lorentz factor of the emitter in the comoving frame is Y' e ~ several, 
the photon interaction angle is at most 2/Y' e . This will reduce the optical depth for pair production. The inferred r (given a same Y) 
should be smaller. For example, for Y' e ~ 2 the optical depth is lower by a factor of ~ 16 and the inferred r decreases by a factor of 
~4. 

A possible pair attenuation exponential cutoff signature may have been observed in the pulse 2 of GRB 060105 with the 
joint Swift-Konus-Wind data (Godet et al. 2007). The observed cutoff is around 600 keV, and the spectral index before the cutoff 
is flat: /3i = 0.67. The observed photon number flux at 1 MeV is 0.5 photons cm -2 s _1 MeV -1 , and the observed duration of 
the pulse is about 20sec. A pseudo redshift z ~ 4 is inferred, which is consistent with the broad pulse profile in the observed 
lightcurves. This burst likely belongs to Case (III) discussed above. If the cutoff feature is real and is indeed due to pair attenuation, 
the requirement that E° th should be smaller than E° c demands Y < 6. Using Eq.dl8b. l |19l >, and &201 , with /3i = 0.67, one can 
estimate r > 10 16 cm, which is consistent with the conclusion drawn from analyzing the Swift X-ray data (Kumar et al. 2007). 
Since the quality of the data is poor, we look forward to the high quality data retrieved by GLAST to finally pin down r in the future. 

We thank the anonymous referee for important remarks and Olivier Godet, Kohta Murase, Enrico Ramirez-Ruiz for helpful 
discussion/comments. This work is supported by NASA under grants NNG06GH62G, NNX07AJ66G and NNX07AJ64G. 
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